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Electron vortex beams carry well-defined orbital angular momentum (OAM) about the propa- 
gation axis. Such beams are thus characterised by chirality features which make them potentially 
useful as probes of magnetic and other chiral materials. An analysis of the inelastic processes in 
which electron vortex beams interact with atoms and which involve OAM exchange is outlined, lead- 
ing to the multipolar selection rules governing this chiral specific electron vortex beam spectroscopy. 
Our results show clearly that the selection rules are dependent on the dynamical state and location 
of the atoms involved. In the most favorable scenario, this form of electron spectroscopy can induce 
magnetic sublevel transitions which are commonly probed using circularly polarized photon beams. 
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Particle vortices, most notably electron vortices (EVs), 
are currently the focus of much interest following the pre- 
diction of their existence on a theoretical basis by Bliokh 
et al. pQ. Such vortices have been created in a num- 
ber of laboratories, using various techniques. 2 -7]. This 
area has recently emerged after much fruitful research 
has been carried out on optical vortices (OVs) over the 
last two decades or so, which has led to a wealth of fun- 
damental knowledge and significant applications [5HTU], 
Both the optical and electron vortices are characterised 
by the singular nature of their wavefronts, with a well 
defined vortex core and quantised orbital angular mo- 
mentum (OAM) about the vortex axis. The general ex- 
pectation is that in all cases the vortex OAM should play 
an important role in the interaction of the vortex with 
matter. However, in the case of an OV, a dipole active 
transition involves exchange of OAM with the centre of 
mass only [TTJ [TJ] , a finding which has been confirmed 
experimentally [HI QI|. The transfer of the OAM of 
an OV to the centre of mass of bulk matter has been 
demonstrated, most clearly in the optical spanner effect 
[Pol Ho] , but as far as optical spectroscopy is concerned 
the development of OAM-based OV beam spectroscopy 
has been hampered by the weakness of optical multipo- 
lar transitions. In contrast, we have recently shown that 
OAM can be transferred efficiently from an EV beam to 
atomic electrons through dipole active transitions [12] [P7] 
and experimentally, a dichroic electron energy loss spec- 
troscopy signal has been demonstrated [3], opening up 
the prospect of chiral specific electron vortex beam spec- 
troscopy (CEVBS) based on OAM selection rules. This 
article seeks to determine the key characteristics of OAM- 
and chiral-related EV beam spectroscopy by analysing 
the OAM selection rules, and specifying the spatial de- 
pendence of the atomic transitions involved. 

The leading interaction between the EV and an atom 
possessing Z electrons is given by the Coulomb interac- 
tion Hamiltonian 
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of the beam electron, the j-th atomic electron and the 
nucleus, all expressed relative to the laboratory frame of 
reference. The transition matrix element between overall 
eigenstates of the atom and vortex system can be written 



as M 



F 



I j, where \I) and \F) are, respec- 
tively, the initial and final unperturbed quantum states of 
the overall system, being products of unperturbed quan- 
tum states of the EV and the the atom: \tpEv) IV'atom)- 
This is justified as the energy of the EV beams inves- 
tigated so far (lOOkV or more) is much larger than the 
energy of the atomic electrons and that of the atomic cen- 
tre of mass. The interesting case in which a low energy 
vortex beam interacts with the atom requires a different 
treatment involving particle exchange, and is not consid- 
ered any further here. In the present case the atomic 
quantum state can be taken as a product of the quan- 
tum state of its nucleus, here taken to also be the cen- 
tre of mass of the atom, and that describing the inter- 
nal electronic state relative to the centre of mass R, i.e. 

IV'atom) = IV'cro(R)) \lpq(ri, ■■■Tj, ...T Z )). 

Without loss of generality, we assume a Bessel EV 
beam of winding number I which is an eigenstate of the 
Schrodinger equation for a free electron in cylindrical po- 
lar coordinates (p vi <p v ,z v ) with the beam axis along the 
z-direction, 
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where kj_ and k z are the transverse and longitudinal com- 
ponents of the wavevector of the vortex beam such that 
k± + = k 2 = ^p^, with E the beam energy, and 
Ji{k±p v ) is the Zth order Bessel function satisfying the 
orthornormalisation condition 



where r„, r„ and R are the position vectors, respectively, 
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Ji(k_ L p v )J l ,(k' J _p v )p v dp v = —5ij>5(kj_ - k'j_), (3) 

fc_L 

which, along with the usual delta function normalisation 
over z gives the normalisation factor of Eq. |2]) . 

CEVBS is concerned with processes in which an inci- 
dent EV mode \k±,l,k z ) lah is scattered by the atom into 
an outgoing EV mode Wj_,l' ,k' z ) iah , with the atom un- 
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dergoing a quantum transition between its internal eigen- 
states. The treatment can be readily extended to more 
general vortex beams since such beams can be repre- 
sented by a linear combination of the Bessel basis modes 
discussed here. As a simplification which is justifiable for 
spectroscopy of condensed matter systems, we shall ini- 
tially assume that the scattering process does not alter 
the state of the atomic centre of mass, which we take as 
stationary. Although this fixed centre of mass assump- 
tion is adequate in the present context of chiral specific 
spectroscopy, the centre of mass coordinate must be re- 
tained as a dynamical variable, as done in the case of 
OVs [TO] HE] [19] when the gross atomic motion is of in- 
terest due to the influence of forces and torques imposed 
by the EV on the atom. 

As in the case of the optical vortex interacting with the 
atom, the transition matrix element for an EV interacting 
with the atom may be reduced to the following form [20] 
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where \i) and |/) are the initial and final states of the 
atom. The effective operator 0> acts on a single elec- 
tron, and emerges in the form 
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where 



F 1 / = Ji(k ±Pv )Mk> ± p v )e i <- l - l '»<'. 
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Note that the first term in Eq. ([lj does not contribute to 
the matrix element by virtue of the orthogonality of the 
initial and final atomic states \i) and |/). 

The chief difficulty in the evaluation of the effec- 
tive operator for the vortex beam-atom interaction in 
Eq. ([5| stems from the fact that the vortex state func- 
tion is conveniently expressed in terms of the laboratory 
frame, while the internal atomic states are customarily 
expressed in spherical coordinates in a frame of reference 
centred on the atomic centre-of-mass of coordinate R. In 
order to overcome this difficulty, the addition theorem of 
Bessel functions [21] [22] can be utilised to represent the 
original EV beam of mode I as a sum of other vortex 
states relative to a shifted frame of reference centred on 
the atomic centre of mass coordinate R. The addition 
theorem reads 
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where a, b and c are three sides of a triangle, and 9 and 
ip the internal angles between sides a and b and b and 
c respectively. Applying this to the triangle formed by 



the position vectors of the vortex, nucleus and atomic 
electron, we identify r c (p c , <f) c , z c ) as the position vector 
describing the vortex electron relative to the centre of 
mass, and after some further algebraic manipulation, we 
find 

oo 

Ji(kxPv) = e- il ^ Ji-p(k ±PR )J p (k ± p c )e^P^e^ 
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(8) 

The EV beam may now be written as the sum of vortex 
modes relative to the centre of mass 
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where the subscripts refer to the reference frame rela- 
tive to which the quantum states are described, and the 
coefficient C;_ p is the overlap integral, given by 



C;-p = Ji- p (kj_p R )e 



i(l-p)4> R 



(10) 



As expected, the above expansion indicates that the 
only vortex mode relative to the centre of mass present 
for an atom located on the beam axis is that for which 
p = I. However, for an atom not situated on the beam 
axis, the strength of the atom-centred vortex modes with 
p = I + 1 and p = I — 1 also become significant when the 
atom is positioned at radial distances of the order of a 
fraction of « 0.1 nm, where a/ i is the first zero of 
the Zth-order Bessel function, i.e. within the first ring of 
the vortex beam. Thus the immediate consequence of 
the shift of the axis is the importance of vortex modes of 
winding numbers different from I, relative to the atomic 
centre of mass frame, as is the case in the context of OV 
beams. This is a manifestation of the extrinsic property 
of the orbital angular momentum of the beam [23 - once 
we transform from the symmetry axis of the vortex to 
one centred on the atomic centre of mass R the OAM 
also changes, as emphasised by Berry [53] 



L z -> L z + e z ■ p R A P, 
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where P is the transverse current in the vortex beam. 
Thus the OAM assumes different values for an off-axis 
atom at different positions relative to the symmetry axis. 

Applying the shifted wavefunctions of Equations ^ 
and (10), the effective operator O 11 (the subscript j will 



henceforth be dropped) relative to the atomic frame can 
be written as follows 
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and 
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with r q (p q ,9 q ,(f> q ) = Tj - R and r c (p c , 6> c , C ) = 
being the internal electronic and vortex beam coordinates 
respectively, both described about the atomic centre, and 



= J P (k± Pc )J p ,{k>^p c )e^-^ 
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We have also isolated the centre-of-mass factor Fr, rela- 
tive to the the atomic frame, from the integral I c involv- 
ing the atomic states |i) and |/). 

In order to evaluate the matrix element in terms of 
multipole excitations, the effective operator must be ex- 
panded in powers of r q . This can be achieved by invoking 
the addition theorem for Bessel functions again in order 
to achieve a separation of the dependence on the atomic 
electron position variable r q from that of the EV r c . This 
is conveniently done by introducing a relative position 
vector s = (r c — r q ). After some algebraic manipulation, 
the integral can be written as follows 
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where F q is given by 

= J p - u (k ±Pq )J p >- u > (fcVJe"^-"'-^')^ , 

(17) 

and 7" is given by 
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The integral J] 1 '" can be evaluated by expressing each 
Bessel function as a coherent superposition of plane 
waves with phase angle dependent on the topological 
charge |21j . The linear momentum transfer wavevector 
is defined as Q(/3) = k/ — k i7 where is the wavevec- 
tor of the plane-wave components of the incident Bessel 
beam, and k/ is that of the outgoing Bessel beam, with 
/3 = 4> — 4>' the relative azimuthal angle between k^ and 
kf. The integral over the vortex beam spatial variables 
may now be evaluated as the Fourier transform of the 
Coulomb potential, leading to the condition u = v!\ 
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We note that 1 /Q 2 (/?) is the familiar kinetic factor arising 
in Coulomb scattering. 

The result for the effective operator of the vortex elec- 
tron beam can now be determined by combining Equa- 
tions (pb, (pi (fl6|, (fl7| and 
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out-of-plane excitations. Eq. (20 1 allows a clear descrip- 
tion of the effect of the EV beam expansion - contained 
in the Fr factors - and its implications for the OAM 
transfer between the EV and the atomic electron - con- 
tained in the F q factors. We illustrate this by considering 
the implications for chiral specific vortex electron beam 
spectroscopy. 

Since the effective operator O 1,1 acts on the electronic 
states of the atoms only through terms containing com- 
ponents of r g , only the terms involving O z and F q are 
relevant. Focusing on O z first, we can write e l ( kz ~ k ^ Zq 
as 1 — i(k z ~ k' z )z q + 0(z 2 ). Thus, the out-of-plane exci- 
tations importantly have no chirality features. 

On the other hand, the term F^~ u p ~ u depends on the 
in-plane components of r q , and contains the phase fac- 
tor e^ p ~ p ^ q which is important for chiral specific spec- 
troscopy. This becomes clear if we consider the simplest 
case of an atom located on the beam axis, in which case 
R = 0. We see that Fr is non-zero only for p = I and 
p' = I', so that the summation over p,p' in Eq. (20) 
amounts only to a single term with p = I and p' = V, and 
Fr = 1. The simplified operator can then be written, in 
ascending powers of p q , 
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In arriving at this result, we have expanded the Bessel 
functions involved in the asymptotic limit of small argu- 
ment as [21] 



Jn(k±p q ) = (-1) n J- n {k±p q ) 



1 / k±p q 
r(n + 1) V 2 



(22) 

where T(z) is the Gamma function. This is equivalent 
to restricting EV spectroscopy to the limit in which the 
transverse momentum of the beam is small compared to 
its axial momentum, a condition often observed in high 
energy electron energy- loss spectroscopy [25] , In such 
cases, the dipole terms in Eq. (21 1 containing the fac- 



tors p q e ±l ^ q Si^±i, when operating on the atomic state, 
causes the magnetic quantum number m of the electronic 
state to increase or reduce by one unit as a result of the 
transfer of one unit of OAM from the EV beam, leading 
to the selection rule l — l' = —m + ml = ±1. This is pre- 
cisely the result obtained by Lloyd et al. [HI [T7] . The 
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FIG. 1. (Colour Online) Transfer of OAM in an atom-vortex 
interaction in the cases when the atom is situated (a) on the 
beam axis, and (b) off axis. In (a) the resulting I' states 
relate directly to an interaction in which I — I' units of OAM 
are exchanged; however it is not possible to determine, for 
example, whether an interaction for which I — l' = 1 is due 
to a dipole interaction, or a higher multipole. In (b), the 
final states I' arise due to a set of transitions in which varying 
quantities of OAM are exchanged with the different p-modes 
the atom "sees". As an example, the p = I — 1 mode may 
be scattered to the p' — p + 2 mode - which may occur in 
a quadrupole or higher multipole interaction - resulting in a 
final state I' — l + l. Thus the final states of the vortex 
beam are not an obvious indicator of the multipolar nature 
of the interaction. The strengths of the modes p and p' are 
dependent on several factors, including transverse coordinate 
Pr and atomic multipolar interaction strengths; thus the bar 
heights are for illustrative purposes only. 



analysis can be extended to higher powers of p q , leading 
to higher multipole excitations involving OAM transfer 
involving multiple units of h. This situation is depicted 
in Fig. Q(a). 

For I — 0, we have A +1 = — A^ 1 ' . Thus, besides the 
phase factor, the effective dichroic operator for a vortex 
beam interacting with an atom is directly comparable to 
the operator associated with the absorption and emission 
of a circularly polarised photon, O ~ (e x ± ie y ) ■ r q — 
x q ± iy q = p q e ±% ^i [I2~| . The dichroic EV beam electron 
energy loss spectroscopy can be seen to be directly equiv- 
alent to the optical circularly polarised dichroic absorp- 
tion spectroscopy by using an I — incident Bessel beam, 
analysing the vorticity of the outgoing beams and com- 
paring the intensity of the exiting V = ±1 beams. Since 
the formal equivalence is at the operator level, the re- 
sult applies to any quantum system with axial symmetry 
about R. In this regard, vortex beam chiral spectroscopy 
is similar to EMCD [26J but would be much more practi- 
cal because only small angle scattering is involved in the 
vortex beam case, so the signal-to-noise ratio should be 
much improved. 

We can now address the physical meaning of the dou- 



tions of this for the chiral-specihc spectroscopy of atoms 
located away from the vortex beam axis. It can be seen 
that the features uncovered above as regards OAM trans- 
fer from the EV to the atomic electron still apply, except 
that now the vortex states with which the atom interacts 
are characterised by the winding number p, not I, and the 
outgoing states after the electronic transition within the 
atom are characterised by p' rather than V . This can be 
understood as a mode broadening effect of the EV beam 
arising from the translation of the reference axis, as de- 
scribed by Eq. ^ and Eq. (10). However, since chiral 
specific spectroscopy is normally conducted with respect 
to the beam axis, summed over atoms at various off-axis 
positions, the spectral changes observed in different OAM 
components of the EV beam can not be exactly related 
to the change in OAM of the atomic electronic system, 
as has been assumed in the case of [3|. The off-axis case 
is illustrated in Fig. 0(b). 

An exception is in the case of an atom whose cen- 
tre of mass is in a pure OAM state, such as found in a 
Bose- Einstein condensate 27, 28J. In that case, the OAM 
states of the atoms would contribute a factor e l ( L ~ L >^ R 
within the matrix element and we must integrate the fac- 
tor Fr with respect to the dynamical variable 4>r. 
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This gives rise to a selection rule for OAM transfer in- 
volving the atomic centre of mass such that Sp — 8L = 61. 
Sp then corresponds to the net OAM change induced in 
the atomic system and so we recover the selection rules 
derived in [T21 E] . One way to understand this result 
is to view the azimuthally delocalized state of the atom 
interacting coherently with the vortex beam. It has been 
shown by Berry that the OAM of the beam is then purely 
intrinsic |24j . If we view the atom interacting with the 
EV beam as a detector, then the conservation of orbital 
angular momentum is dependent on the nature of the 
detector, and whether it has the cylindrical symmetry of 
the vortex beam. This argument could lead to a general 
test as to whether or not a system is coupled coherently 
to the EV. 

In summary, we have carried out a comprehensive anal- 
ysis of OAM transfer in inelastic atom-vortex interac- 
tions, and derived the effective operator exhibiting quan- 
tised OAM transfer via multipole electronic excitation of 
the atom. When the atom is on-axis, spectroscopy with 
the I = beam is shown to be equivalent to magnetic 
circular dichroism, however, for the cases when the atom 
is off-axis, the extrinsic nature of the OAM of the vor- 
tex beam leads to a complicated situation in which the 
change in OAM of the vortex beam does not directly 
correlate with the change in OAM of the atom. As such, 
simple application of global OAM conservation is not the 
correct approach in analysing such interactions, and the 
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implication for chiral spectroscopy is that extreme care 
must be take in experimental analysis, with the relative 
strengths for the multipole interactions taken into ac- 
count. The dependance of the interaction on coherent 
coupling with the atomic state is noted, and this suggests 
a general test for coherence of the vortex interactions. 
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